QUANTUM PROBABILITY 
APPLIED TO THE DAMPED HARMONIC OSCILLATOR 
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Abstract. In this introductory course we sketch the framework of quantum proba- 
bility in order to discuss open quantum systems, in particular the damped harmonic 
oscillator. 
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1. The framework of quantum probability 

Noncommutative probability theory (or 'quantum probability') generalises Kol- 
mogorov's classical probability theory in a way that allows the inclusion of quantum 
mechanical models. For a discussion of its motivation we refer to |KiiM| in this se- 
ries. Basic sources on quantum probability outside the present series are |Biaj . |Davj . 
Gudj . [Mac], |Mey| , |Neuj . [Par], |Var j . An independent introduction is given 



Hoi 

here 

Making probability noncommutative. In the last two decades a succesful strat- 
egy has become popular in mathematics: the generalisation of classical mathematical 
structures by noncommutative algebraic constructions. The most widely known ex- 
ample where this strategy was applied is doubtlessly the noncommutative version of 
geometry, as explained in the imaginative book of A. Connes ( |Conj ). There the clas- 
sical structures of a topological space and of a differential manifold are the pillars on 
which the K-theory of C*-algebras and a variety of coho mo logical algebras are built. 
Another application is the field of 'quantum groups', where the classical structure of 
a Lie group leads into new areas in the theory of Hopf algebras. 

However, the oldest case by far is von Neumann's and Segal's 'noncommutative 
integration theory', which has developed into noncommutative measure theory and 
probability theory. 

The general strategy consists of the following three steps. 

l 
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(1) Encode the information contained in the classical structure into an appropri- 
ate algebra of functions on it. 

(2) Characterise the resulting algebra axiomatically. One of the axioms will be 
commutativity. 

(3) Drop the commutativity axiom. 

Classical probability. Let us apply this strategy to the structure of a probability space. 

We remind the reader that a probability space is a triple (Q, E, P), where Q is a set, 
E is a a-algebra of subsets of Q, containing Q itself, and P is a u-additive function 
E -> [0, 1] with the property that P(fi) = 1 (|Kol)). 

In applications Q is interpreted as the set of all possible outcomes of a certain 
stochastic experiment. E consists of 'events', statements about the outcome of the 
experiment that can be tested by observation. When E is such an event, then F(E) 
is the probability that E will occur. 

Applying the strategy. Step 1. We choose to consider the algebra L°°(Q, S,P), con- 
sisting of all bounded measurable functions /: Q — > C, where two such functions 
/ and g are identified if / — g vanishes P-almost everywhere. On this algebra we 
consider the linear functional p given by 



We have chosen measurable functions because we want to encode a measurable struc- 
ture. (Had we been interested in the topological structure, we would have chosen an 
algebra of continuous functions.) 

Bounded functions are an appropriate choice since they allow for unlimited multi- 
plication of elements, while keeping (p well-defined. 

The identification of functions which are almost everywhere equal is a technical 
simplification, standard in integration theory. 

Now we must check whether all the relevant information in (Q, E,P) has been faith- 
fully encoded. Clearly, the triple (fi, E, P) determines L°°(Q, E, P) uniquely. In the 
converse direction, we recover a cr-algebra E by putting 



which, however, is not isomorphic to E, since we have identified functions that are 
equal almost everywhere. In fact, E is the measure algebra, that is the quotient of E 
by the equivalence 




E := { p E L°°(tt, E, P) | p = p 2 = p* } , 



S ~T, meaning: F((S \ T) U (T \ S)) = . 



This simplification is a gain rather than a loss. 

Finally the probability measure P is regained by putting 



P : E -> [0, 1] : p ^ <p(p) . 



Step 2. L°°(Q, E,P) is characterised as a commutative von Neumann algebra. A 
few definitions should now be given. 
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Let Ti be a Hilbert space and let Ai, A 2 , A 3 , • • • be a sequence of bounded operators 
on H. This sequence is said to converge to a bounded operator A in the strong 
operator topology if for all tp G H: 

lim || A n %l) - Atp || = . 

n^oo 

It increases to A if, moreover, Aj < Aj +1 in the sense that Aj +1 — Aj is a positive 
operator. 

A von Neumann algebra A is an algebra of bounded operators on some Hilbert 
space Ti which is closed in the strong operator topology. We shall always assume 
that A contains the identity operator 1, and we only consider separable Hilbert 
spaces. 

A state on a von Neumann algebra A is a linear functional ip : A — > C mapping 1 
to the number 1, and which is positive in the sense that p(A*A) > for all A G A. If 
ip(A*A) = only for A = 0, then <p is called faithful. If \im\{p(A\) = p(A) for every 
net {A\) of positive operators increasing to A, then p is called normal. 

Finally, for a bounded function / on a probability space (Q, S, P) we denote by Mj 
the operator of multiplication by / on the Hilbert space L 2 (Cl, S, P). 

Proposition 1.1. Let E, P) 6e a probability space. Then the algebra 

A:= {M f \ f G L°°(fi,E,P)} 

is a (commutative) von Neumann algebra of operators on L 2 (fl, E,P), and i/ie map 
<£> : Mf i— > J /dP «5 a faithful normal state on A. Conversely, every commutative 
won Neumann algebra A with a faithful normal state p is of the above form for some 
classical probability space. 

Proof. We only prove the first part of the theorem. The point is to show that A is 
strongly closed. So let (fx) be a net of L°°-functions such that Mf x tends strongly to 
some bounded operator X on Ti := L 2 (Q, E, P), that is for all ip G Ti. we have 

L 2 -hmM> = A^. 

A 

Without loss of generality we may assume that || X \\ = 1. We must show that 
X = M f for some / G L°°. Put / := XI. Then for all g G L°° we have 

Xg = L 2 -\imf x g = L 2 -\imM 9 f x = M g (L 2 -lim/ A • l) = M g Xl = M g f = fg . 

A A \ A / 

Now let the event E £ for e > be defined by 

E £ :={iuen \ \f(u)\ 2 >l + s} . 

Then, since || X || < 1, 

F(E £ ) = || l Ee f > || X1 Ee f = II fl Ee || 2 = / |/fdP > (1 + £ )P(£ £ ) , 

and it follows that F(E £ ) = 0. Since this holds for all e > 0, we have |/| < 1 almost 
everywhere with respect to P. So / G L°°(Q, S,P). Finally, since the operators X 
and Mf are both bounded and coincide on the dense subspace L°° of H, they are 
equal. □ 
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Step 3. We now drop the commutativity requirement to arrive at the following 
definition. 

Definition. By a noncommutative probability space we mean a pair (A, ip), where A 
is a von Neumann algebra of operators on some Hilbert space 7i, and ip is a normal 
state on A. If ip is faithful, the probability space is called non- degenerate. 

Events and random variables. Let us carry some important concepts of proba- 
bility theory over from the structure (Q, E, P) to the generalised probability space 
(A,<p). (Comparable discussions are found in [Mac], |KiiM| . (Par], |Mey| .) 

Classically an event is an element S of E. In Step 1 of the preceeding section this 
is replaced by the projection l s in L°°(fl, E, P). In Step 3 the concept of an event 
is generalised to that of an arbitrary orthogonal projection in A, that is an element 
satisfying E 2 = E = E*. The state p associates to this event the probability <p(E). 
The operator is the impossible event, and 1 is the sure event. Two events E and 
F are called compatible if EF is also an event, equivalently E and F commute: 

EF = (EF)* EF(— (EF)*) = FE => (EF) 2 = E 2 F 2 = EF . 

If this is the case, then the event EF stands for the occurrence of both E and F and 
the event E V F := E + F — EF for the occurrence of either E or F or both. If 
EF = 0, then the occurrences of E and F exclude each other. So mutually exclusive 
events are described by orthogonal subspaces of Tt. 

A classical random variable on a probability space (Q, E, P) is a measurable func- 
tion X from Q to some other measure space (fi',E'). Such a function induces an 
embedding of E' into E given by 

J x : S i-h. X-\3), 

containing the same information as X itself. The probability distribution of X is 
given by 

F x := P o J x : S i-> P(X~ 1 (S)) . 

Let us see what our program does with this structure. In Step 1 the embedding Jx 
is replaced by the mapping 

jx : L°°{a, E\ F x ) - L°°(Q, E, P) : / h- / o X , 

the natural extension of the map 1$ i— > E(S) := lx- 1 (s) to an (injective) *-homomorph- 
ism L°°(n', E',P X ) -> L°°(Q, E,P). Still the projection E(S) stands for the classical 
event X^ X (S) that the random variable X takes a value in S G E'. 
In Step 3 this is now generalised to the following notion. 

Definition. By a generalised random variable on a noncommutative probability space 
(A, (p) we mean a *-homomorphism from some other von Neumann algebra B into A 
mapping lg to 1^. The probability distribution of j is the state ■?/> := ip o j on B. 

We denote this state of affairs briefly by 

j:(B^)^(AM- 
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If B is commutative, say (B,ip) = L°°(Q', E',P'), then the random variable j is said 
to take values in Q', and j can be written 

3(f) = I fWE(dX) , 
Jn> 

where E denotes the projection-valued measure given by 

E(S):=j(l s ), (SEE'). 

In the particular case that Q' = R, j determines a unique self-adjoint operator on 
the representation space H of A: 

Theorem 1.2 (Spectral Theorem, von Neumann). There is a one-to-one correspon- 
dence between self-adjoint operators A on a Hilbert space TC and projection-valued 
measures E : E(R) — ► B(fH) such that 

A = [ \E(d\) . 
Jr 

When E(S) G A for all S in the Borel a-algebra E(R), then A is said to be affiliated 
to A. Moreover: 

Theorem 1.3 (Stone's Theorem). There is a one-to-one correspondence between 
strongly continuous unitary representations 1 1— > U t of the abelian group R into A and 
self-adjoint operators A affiliated to A such that 

U t = e itA . 

Here the right hand side is to be read as the strongly convergent integral 

e itA ._ f e it\ E ( dX ) } 



where E is given by the spectral theorem. If we put e t : R — > C : x i— > e ttx , then the 
connection with j can be written 

3{e t ) = e ltA . 

So altogether we can characterise a real-valued random variable or observable in any 
one of four ways: 

(1) by a self-adjoint operator A affiliated to A; 

(2) by a projection- valued measure E in A; 

(3) by a normal injective *-homomorphism j : L°°(R, E(R),P) — > A; and 

(4) by a one-parameter unitary group (U t )tm in A. 

Interpretation of quantum probability. It is a surprising fact that nature — at 
least on small scales — appears to be governed by noncommutative probability. 

Quantum probability describes manipulations performed on physical systems by 
certain mappings between generalised probability spaces called operations. (The same 
has been said about classical probability see, e.g. |Kamj .) These mappings will be 
treated in some detail in Section El The generalised random variable which we just 
saw is such a mapping. It represents the operation of restricting attention to a 
subsystem. Another such mapping is the conditional expectation, describing the 
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immersion of a physical system into a larger one. Yet other operations are the time 
evolution and the transition operator: they represent the act of waiting for some time 
while the system evolves on its own, or in interaction with something else respectively 
At the end of a chain of operations we land in some probability space (A, <p), and 
we need a way to interpret it in terms of the outcome of the physical experiment 
which is being described. The rules for interpretation are as follows. 

• Some of the orthogonal projections in A have an interpretation as observable 
events. 

In the language of Mackey (JVlacj) events may be considered as questions, which 
can be asked to the system. It answers by saying 'yes' or 'no'. 

• The experiment can be repeated arbitrarily often. Each time we are free to 
choose new questions, that is we are allowed to adjust observation equip- 
ment without our whole experiment becoming a different one. (This is the 
main distinction with the stochastic experiments envisaged by Kolmogorov in 

• Compatible questions can be asked together in the same trial. 

• Incompatible questions can be asked in different trials. 

• Inside one single trial it is sometimes possible to ask incompatible questions 
one after the other. The order will then influence the probabilities: if the ques- 
tions Ei, E 2 , . . . , Eft are asked in each trial and in this order, the asymptotic 
fraction of the trials in which they are all answered 'yes' is 

(p(EiE 2 ■ ■ ■ E k _iE k E k _i ■ ■ ■ E 2 Ei) . 

Remark. It is sometimes difficult to say where the operation ends and the observa- 
tion begins. For example, posing the question E after the question F can alternatively 
be viewed as an operation Z°°({0, 1} x {0, 1}) — > A, followed by the observation of 
the compatible events {1} x {0, 1} and {0, 1} x {1}. 

The quantum coin toss: 'spin'. The simplest noncommutative von Neumann 
algebra is M 2 , the algebra of all 2 x 2 matrices with complex entries. And the simplest 
noncommutative probability space is (M 2 , |tr), the 'fair quantum coin toss'. 

The events in this probability space are the orthogonal projections in M 2 : the 
complex 2x2 matrices E satisfying 

E 2 = E = E* . 

Let us see what these projections look like. Since E is self-adjoint, it must have two 
real eigenvalues, and since E 2 = E these must both be or 1. So we have three 
possibilities. 

(0) Both are 0; that is E — 0, the impossible event. 

(1) One of them is and the other is 1. 

(2) Both are 1; that is E = 1, the sure event. 

In Case (1), E is a one-dimensional projection satisfying 



tr £" = + 1 = 1 and det E = • 1 = 



As E* = E and tr E = 1 we may write 
E = 



1 + z x — iy 
x + iy 1 — z 



with (x, y, z) G 



Then det E = implies that 

H(l-z 2 )-(x 2 +y 2 )) 







2 , 2 , 2 

x + y + z 



1 . 



So the one- dimensional projections in M 2 are parametrised by the unit sphere 5*2. 
Notation. For a = (ai, a 2 , 03) G M 3 let us write 



a a 



03 a x — za 2 
a x + ia 2 —03 



ai<Ti + a 2 cr 2 + a 3 a 3 



where &\ , o<i and (T3 are the Pauli matrices 





"0 1" 




"0 -i" 




"1 " 


0i := 


1 


, CT2 : = 


z 


, 03 : = 


-1 



We note that for all a, b G M 3 we have 

cr(a)a"(6) = (a, 6) • 1 + ia(a x 6) . 

Let us write 

£(a):=±(l + a(a)), (||o|| = l). 
In the same way the possible states on M 2 can be calculated. We find that 



(f{A) = tr (pA) where p = p(a) 



(j(a)), 



< 1 . 



The situation is summarised by the following proposition. 



(1.1) 
(1.2) 

(1.3) 

as in 



Proposition 1.4. The states on M 2 are parameterised by the unit ball in '. 
fjl.3|) . and the one- dimensional projections in M 2 are parametrised by the unit sphere 
as in (|1.2|) . The probability of the event E(a) in the state p(b) is given by 

tr (p(b)E(a)) = 1(1 + (a,b)) . 

The events E(a) and E(b) are compatible if and only if a = ±6. Moreover we have 
for all a G S 2 : 

E(a) + E(-a) = 1 , E(a)E(-a) = . 

Interpretation. The probability distribution of the quantum coin toss or l qubiV is 
given by a unit vector b in three-dimensions. For every a on the unit sphere we can 
say with probability one that of the two events E(a) and E(—a) exactly one will 
occur, E(a) having probability |(1 + (a, b)). We therefore have, for each direction a, 
a classical coin toss with probability for heads equal to |(1 + (a,b)). The coin tosses 
in different directions are incompatible. 

Particular case: the quantum fair coin is modelled by (M 2 , |tr). 

The quantum coin toss is realised in nature: the spin direction of a particle with 
total spin \h~ behaves in this way. 



Positive definite kernels. In this section we introduce a useful tool for the con- 
struction of Hilbert spaces, used heavily in quantum probability. 

Let S be a set and let K be a kernel on S, that is a function S x S — > C. Then X 
is called positive definite if for all n G N and all n-tuples (Ai, . . . , A n ) G C n we have 

n n 

t=i j=i 

Theorem 1.5 (Kolmogorov's dilation theorem). Let K be a positive definite kernel 
on a set S. Then up to unitary equivalence there exists a unique Hilbert space TC and 
a unique embedding V : S — > TC such that 

V x , veS : (V(x),V{y)) = K{x,y) (1.4) 



and \/V(S) =W . (1.5) 

A map V : 5 — ► TC is called a (Kolmogorov) dilation if ()1.4j) holds. It is called 
minimal if (ll.5j) holds. 

Proof. Consider the space £ of all functions 5 — > C with finite support. Then £ 
becomes a pre-Hilbert space if we define the (pre-)inner product 

(A,/i) := ^2^2\(x)K(x,y)fM(y) . 

x&S ydS 

Dividing out the null space 

AT := { A e £ | (A, A) = } 
and forming the completion TCk of £/Af, we let Vk '■ S — > TCk be given by 

V K {x) :=5 X +Af. 

Then for all x,y G S: 

(V K (x), V K {y)) = (5 X + Af,5 y + M) C /m = ($x, S y )c = K(x, y) . 

Now let V : S — ► TC be a second minimal Kolmogorov dilation of K. Then we define 
a map 

U :£^TC: H x ) v ( x ) . 

This map vanishes on A^: for A G Af we have 

2 

|2 " 



C/ A 



X)A(x)y(x) 



= E E j/)A(y) = (A, A> £ = . 

So Z/o may be considered as a map C/Af — > 7Y. By the same calculation we find that 
£/o is isometric. Since \/ V(«S) is dense in TC and V Vr-(«S) is dense in TCk, U extends 
to a unitary map U : TCk TC mapping Vk{x) to V(x). □ 

Examples 1.6. 

(a) Let S be any set and let K(x,y) := 8 XjV . Then 7i = l 2 (S) and V maps the 
elements of S to the standard orthonormal basis of TC. 
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(b) Let S := Hi x H 2 , the Cartesian product of two Hilbert spaces Hi and H 2 . 
Let 

K((ipi,ijj 2 ), (Xi,X2)) := (0i, Xi) • (02,X2> • 
Then H = Hi®H 2 , the tensor product of Hi and H 2 , and V(ipi, i> 2 ) = ^i®^2- 

(c) Let S be a Hilbert space; call it K, for the occasion. Let K(ip,x) '■= (0, x) 2 - 
Then H is the symmetric tensor product /C © s /C and V(ip) = ip © 0. 

(d) Let 5 := /C as in Example (c). Let K be the positive definite kernel 

Then the Kolmogorov dilation is the Fock space TiJC) over IC, defined as 

F{K) :=C©/C©±(/C© S /C) ©±(/C© s /C© s /C) © 

and V{ip) is the so-called exponential vector or coherent vector 
Exp (0) := 1 © © (0 © 0) © (0 © © 0) © • • • 

(e) Let 5 = R and let K : R x R -> C be given by 

i) := e -il«-t|+M«-*) ? > o, u e M). 

The Kolmogorov dilation of this kernel can be cast in the form 

{ (r)-iui)(x-t) if T < /■ 
if X > t. 

2. Some quantum mechanics 

Quantum mechanics is a physical theory that fits in the framework of noncom- 
mutative probability, but which has much more structure. It deals with particles 
and fields, using observables like position, momentum, angular momentum, energy, 
charge, spin, isospin, etc. All these observables develop in time according to a certain 
dynamical rule, namely the Schrodinger equation. 

In this section we shall pick out a few elements of this theory that are of partic- 
ular interest to our main example: the damped harmonic oscillator considered as a 
quantum Markov chain. 

Position and momentum. Let us start with a simple example: a particle on a 
line. This particle must have a position observable, a projection valued measure on 
the Borel a-algebra S(R) of the real line R: 

E : S(R) -> B{H) . 

The easiest choice (valid when the particle is alone in the world and has no further 
degrees of freedom) is 

H := L 2 (R) ; 
E(S) : ip i-> l s - i) . 

In this example the Hilbert space H naturally carries a second real-valued random 
variable in the form of the group (T t ) t6 R of spatial translations: 

{T t ^){x) :=^x-M) , (2.1) 
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according to the remark following Stone's theorem ( Theorem 11.3)1 . This second ob- 
servable is called the momentum of the particle. The constant h is determined by 
the units of length and of momentum which we choose to apply. The associated 
self-adjoint operators are Q and P given by 

{Qiji){x) = xip(x) ; 

(P^)(x) = -ih^(x) . (2.2) 

Just as we have T t = e~ ltp , it is natural to introduce S s := e ts ® whose action on 
H is 

S s i/j(x) := e isx ijj(x) . (2.3) 
The operators P and Q satisfy Heisenberg's canonical commutation relation (CCR) 

[P, Q] = -ih-l. (2.4) 
A pair of self-adjoint operators (P, Q) satisfying (|2.4jl is called a canonical pair. 

Representations of the canonical commutation relations. What kinds of canonical 
pairs are there? 

Before this question can be answered, it has to be reformulated. Relation (|2.4j) is 
not satisfactory as a definition of a canonical pair since the domains on the left and 
on the right are not the same. Worse than that, quite pathological examples can be 
constructed, even if (|2.4jl is postulated to hold on a dense stable domain, with the 
property that P and Q only admit unique self-adjoint extensions ( jReSj ). 

In order to circumvent such domain complications, Weyl proposed to replace (|2.4J) 
by a relation between the associated unitary groups (T t ) and (S s ), namely: 

T t S s = e~ ihst S s T t , (s,teR). (2.5) 

It was von Neumann's idea to combine the two into a two-parameter family 

W(t,s):=ef st T t S s , (2.6) 

forming a 'twisted' representation of M 2 , as expressed by the Weyl relation: for all 
s, t,u,v G M, 

W(t, s)W(u, v) = e-% {tv ~ su) W(t + u,s + v) . (2.7) 

This relation captures the group property of T t and S s together with the relation 
(|2~Hjl . Formally, 

W(t, s) = e i{sQ - tp) . 

We shall call the representation on L 2 (R) of the CCR given by (J2H), (Q, O and 
()2.6|) the standard representation of the CCR. 

Here and in the rest of the text we shall follow the quantum probabilist's convention 
( |Mey| ) , namely that 

h = 2 . 

Theorem 2.1 (von Neumann's Uniqueness Theorem). Let (W(t, s)) t R be a strongly 
continuous family of unitary operators on some Hilbert space 7i satisfying the Weyl 
relation (|2.7jl . Then 7i is unitarily equivalent with L 2 (IR) ® K,, such that W(t,s) 
corresponds to Ws(t,s) (g> 1, where Ws is the standard representation of the CCR. 
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Proof. Let W : R 2 — > W(7Y) satisfy the Weyl relation (|2.7|) . For each integrable 
function / : R 2 — > C with J J s)|d£ ds < oo, define a bounded operator on 
TC by the strong sense integral 



/OO /»00 
/ /(*,sW*,s)d*<fc. 
-oo J — oo 



We find the following calculating rules for such operators A(f) and their kernels f: 
A(f)+A(g) = A(f + g); 



A(f)* = A(f), where f(t,s):=f(-t,-s) ; 
A(f)A(g) = A(J*g). 

Here the 'twisted convolution product' * is defined by 

POO POO 

(f*g)(t,s):= / / e-^-^f(t-u,s-v)g(u,v)dudv. 



oo J — oo 



Moreover we claim that an operator (on a nontrivial Hilbert space) can have at most 
one kernel: 

A(f) = =>► H = {0}oif = 0. (2.8) 
Indeed, if A(f) = then we have for all a, b G M, 

/*OC /»oo 

= W(a, 6)M(/) W (a, 6) = / / e 2l{as - bt) f(t, s)W(t, s)dt ds . 



oo J — oo 



Applying the linear functional A t— > (^>, A0) with </?, ^ e to both sides of this 
equation, we find that for all <p,i/; E TC the (integrable) function 

(M)^/(M)(^iy(M)^> 

has Fourier transform 0. By the separability of TC, either W(t, s) — for some (£, s), 
(that is W = {0}), or f(t, s) — for almost all (t, s). 
The key to the proof of uniqueness is the operator 



1 /"OO />00 

E := — / e^+^V^dtds 

" J —oo J —oo 



It has the remarkable property that for all a, b G R, £W(i, s)i? is a scalar multiple 
of£: 

EW{a,b)E = e-^ a2+b2 )E . (2.9) 
Indeed, has kernel g(t,s) := ^e~2(* +s ) ; and the product W(a,b)E has kernel 

h(u,v) := i e -^-M . e -K( a -") 2 +( 6 -*) 2 ) . 

7T 
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So EW (a, b)E has kernel 

( g *h){t,s) = [ [ e- i{tv - su) g(t-u,s -v)h(u,v)dudv 

e -i(tv-su) e -^((t-u) 2 Ms-v) 2 ) e -i(av-bu) e -^((a-u) 2 Hb-v) 2 ) du ^ 
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7T 2 

J_ e -i(a 2 + b2 ) e -!(* 2 + s2 ) / / e -(«-K*+^+«+*))"-(«-i(*+«+«+^))^ u ^ 
7T 2 





-e 

7T 



which proves 1)2.9}) . 

We conclude that E* = E (since g = g), E 2 = E (putting a = 6 = in 1)2.9}) ) . 
and that EAE = CE, where .4. is the von Neumann algebra generated by the Weyl 
operators. So .E is a minimal projection in A. Denote its range by K,. Then we have 
for all tp, ip G K and all t, s,u,v G K: 

= (<p,EW{-t,-s)W{u,v)Eilj) 
= e l(tv - su) (<p } EW(u-t,v-s)E^} 
= ei{ tv- a u) e -^((u-tr+(v-sr) ^ E ^ 

Therefore the map 

V : R 2 x K -> H : ((t, s), h-> e^ {t2+s2) ^(t, s)p 
is a Kolmogorov dilation (cf. Section QJ) of the positive definite kernel 

K : (M 2 x /C) x (R 2 x /C) -> C, ((t, s), (ix, v), if>) ^ e {t - ls)(u+iv) (<p, iff) . (2.10) 

By explicit calculation you will find that E$ is the orthogonal projection onto the 
one-dimensional subspace spanned by the unit vector Q(x) := 11/7(2), where 

jlx) := -^=e 2 . 



So in the standard case the dilation is 

V s : R 2 -> L 2 (R) : (t, s) i-> e^ +s%) e~i ts e isx Vt{x - 2t) . 

By Kolmogorov's Dilation Theorem, there exists a unitary equivalence U : L 2 ( 
/C — > 7^ such that for all a,fteK and ip G /C: 

C/(W<j(a, 6)fi ® V>) = W(a, b)ip , 

and therefore for all a, 6 G R: 

W(a,6) = C/(W s (a,6) ® l)*/" 1 , 
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provided that the range of V is dense in Ti.. Let C denote the orthogonal complement 
of this range. Then L is invariant for the Weyl operators; let Wo(t, s) be the restriction 
of W(t, s) to C. Construct E := A (g) in terms of W in the same way as E was 
constructed from W . Then clearly E < E, but also E J_ E. So E = A (g) = 
and by (jZHJ) we have £ = {0}. □ 

Exercise. Calculate the minimal projection Eg in the standard representation. 

Energy and time evolution. The evolution in time of a closed quantum system 
is given by a pointwise strongly continuous one-parameter group (a t )teR of ^auto- 
morphisms of the observable algebra A. 

Like in the case of a particle on a line, for a finite number n of (distinguishable) 
particles in <i-dimensional space we take A = B(7i) with Ti = L 2 (W ad ). Since all 
automorphisms of this algebra are implemented by unitary transformations of Ti, the 
group (a t ) is of the form 

a t {A) = UtAU^ 1 . 

It is possible to choose the unitaries so that 1 1— > U t is a strongly continuous unitary 
representation R — > IA(H). We denote its Stone generator by H/h: 

U t = e ltH / h . 

The self-adjoint operator H corresponds to an observable of the system of parti- 
cles, called its energy. The operator H itself is known as the Hamilton operator or 
Hamiltonian of the system. As the Hamiltonian commutes with the time evolution 
operators, energy is a conserved quantity: 

a t (H) = U t HUt l = H . 

The nature of a physical system is characterised by its dynamical law (a term of 
Hermann Weyl, see |Wey| ) . This is an equation which expresses the Hamiltonian 
in terms of other observables. For n interacting particles in R d in the absence of 
magnetic fields the dynamical law takes the form 



nd . 

H = J2^—- p ! + v (QuQ2,...,Q 



i=i 



nd I 



Hi) 



for some function V : (R d ) n — > R, called the potential. The positive constants rrik, k = 
1, • ■ ■ , n are the masses of the particles. (Incidentally we put k(j) := 1 + [(J — l)/d], 
where [ ] denotes integer part, in order to attach the same mass to the coordinates 
of the same particle.) 

Free particles. If V — 0, then U t factorises into a tensor product of nd one-dimensional 
evolution operators, all of the form 

U t = e UH l h = e^ p2 . 

Since the Hamiltonian H = P 2 /2m now commutes with P, momentum is conserved: 

a t (P)=P. 
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On a formal level the time development of the operator Q is found by solving the 
differential equation 

^M^jfiQUr^^MQ)}, (2.11) 

a solution of which is 

(*t(Q) = Q + -P. 

m 

According to the Uniqueness Theorem the canonical pairs (P, Q) and (P, Q + —P) 
are indeed unitarily equivalent. So we expect the evolution of the Weyl operators to 
be the following: 

ixP+iyQ\ _ -ixP+iy(Q+j-P) 



a t (W(x,y))=a t ( e ~- p +^) 



e -i( X -l-y)P + iyQ = W [ x _ Ly 



m 



Proposition 2.2. Let P := —ih^ denote the momentum operator on 7i := L 2 
and let W : 1R 2 -> U(H) be given by Let 

TT i-^P 2 



Then 



U t W(x,y)U t - 1 = W(x--y,y 

m 



Proof. From the definitions of T t and Eg it follows that for all measurable sets Bel 
and all tGl: 

T t E Q (B)Tf l = E Q {B + ht) . 

By the uniqueness theorem irreducible representations of the CCR have the symmetry 
Q — > P, P — > — Q. So we also have the exchanged imprimitivity relation 

V fle ECR)V v6R : S y E P (B)S y 1 = E P (B + hy) . 

Hence for all y, t G R, 

SyU^S- 1 = S y (y°° e-^ x2 E P (d\)^ Sy 1 

oo 

e -^M~ h yf E P (d\) 

-oo 

th .,2 

= 17-t • T_± v ■ e~^ v . 
Multiplying by C/j on the left and by S y on the right we find 

U t W{^y)U- 1 = UtSyU- 1 = e~ l ^y 2 T_j_ y S y = W {-Uj ,yj . 

As T x commutes with U t we may freely add (x, 0) to the argument of W, and the 
proposition is proved. □ 
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By imposing some state tp on A = B(L 2 (R)), all stochastic information on the 
model (A, <p, at) can be obtained from the evolution equation a t (Q) = Q + ^P. For 
example, at large times t the random variable \a t (Q) approaches ^P in distribution, 
provided that ip does not favour large Q values too much. So a position measurement 
at a late time can serve as a measurement of momentum at time 0. This puts into 
perspective the well-known uncertainty principle for position and momentum at equal 
times. 

The Schrddinger picture and the Schrddinger equation. The type of description of a 
system given so far, namely with random variables moving in time, and the state 
ip given once and for all, is called the Heisenberg picture of quantum mechanics. 
In probability theory this is common usage, and we shall adopt it also in quantum 
probability. 

However, quantum mechanics is often thought of in a different way, where one lets 
the state move, and keeps the operators fixed. This is close to Schrddinger's 'wave 
mechanics', and is therefore called the Schrddinger picture: 

If we take for p> a pure (that is extremal) state on the algebra A = B(H) where H 
is, say, L 2 (R nd ): 

ip(A) = (4>,A4>) , (^en,n\\ = i), 

then we can express all probabilities at later times t in terms of the wave function 

tp(x 1 , . . .,x nd ;t) := ([/"VX^i, • • -,Xnd) ■ 

This wave function satisfies the Schrddinger equation, a partial differential equation 
reflecting the dynamical law: 

-ih-^ijj(x 1 ,...,x nd ;t) 

nd j q 2 

= ~9 h2 aT^l' • • • ' X ^ t) + V(x 1: . . . , x nd )ip(x 1: . . . , x nd ; t) . 

zm k (j)n aXj 

If E is an orthogonal projection in 7i, then the probability of the associated event 
can be calculated in the Schrddinger picture by 

<p(a t (E)) = (^,U t EU t - 1 ^) = (Uf^EUr 1 ^) 

= / ip t (xi, ■ ■ ■ ,x nd )(Eip t )(x 1 , . . . ,x nd )dx!, . . . ,dx nd . 

The harmonic oscillator. A harmonic oscillator is a canonical pair (Q, P) of ob- 
servables that under time evolution (a t )tm performs a rotation such as 

a t{Q) — Q cos t + P sin t ; 
a t (P) = —Q sin t + P cos t . 

Since rotation in the plane is symplectic (preserves the area two-form), this evolution 
respects the canonical commutation relation QP — PQ = ih-1. So by the Uniqueness 
Theorem it determines (up to a time-dependent phase) a group of unitary transfor- 
mations (U t )teR of the Hilbert space on which it is represented. (For example, Uk is a 
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unitary transformation of L 2 (R) that sends Q into P and P into — Q in the standard 
representation (|2.2jl . This is the Fourier transform.) 

Making a formal calculation as in ()2.11|) by differentiating the equality 

a t (A) = jtm Ae -itB/n 

we find that a Hamiltonian of the form 

H = \{P 2 + Q 2 ) (2.12) 
can be expected to generate such a rotating evolution. 

The textbook treatment of the harmonic oscillator (e.g. [Han]), follows the elegant 
algebraic reasoning of Dirac, who rewrote the Hamiltonian (j2.12j) as 

H = \{Q - iP){Q + iP) + \i[P, Q] =: ha*a + \% ■ 1 . 

The operators a and a* are then seen to lower and raise the eigenvalue of H, and are 
called the annihilation and creation operators. 

Here we choose to proceed more analytically, seizing the opportunity to introduce 
techniques which will be useful again later on for the treatment of free quantum fields 
and the damped oscillator. 

Our goal is to describe H and U t explicitly. 

Heisenberg's matrix representation. First we note that, since ct t has period 2tt, the 
differences between spectral points of H must be multiples of h. On the grounds of 
fl2.12j) we suspect that H is bounded from below, so let us try 

sp(H) =fm + c. 

We take as our Hilbert space Tin '■= I 2 (N, ^r) with the Hamiltonian given by 

(m)(n) = (hn + c)$(n) . 

The subscript 'if' indicates that on this space we wish to stage matrix mechanics of 
the Heisenberg type. If we define on Hh the 'product' or 'coherent' vectors 

tt(z) := (l,z,z 2 ,z 3 ,---), (zeC), 

then our intended time evolution takes the form 

U t H 7r(z) = e itc/h Ti (e tt z) . (2.13) 

Now we want to represent a canonical pair (P, Q) in this space, or equivalently, Weyl 
operators W(z), that rotate in the same way: U t W(z)U^ 1 = W(e zt z). We note that 

n=Q 

so that we have here another dilation of the positive definite kernel ()2.10j) used in the 
proof of the Uniqueness Theorem. An irredicible representation of the CCR is close 
at hand. Put: 

W h (z)it(u) = e~~ zu ^ A \{u + z), (z,ueC) . 
These operators satisfy the Weyl relation 

W H (w)W H (z) = e- ilm{Wz) W H (w + z) , (2.14) 
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the same as ()2.7|) if we identify W(t, s) with Wuif+is). Clearly we have also obtained 



U t W{z)Ut X = W (e u z) 



(2.15) 



Let us summarise, again replacing h~ by 2. 



Proposition 2.3. The Heisenberg representation of the Harmonic oscillator is given 
by 



l 2 (N, 



1 



7?! 



u t H *(z) 



Hh 

(H H $)(n) = (2n + l)0(n) 
W h {z)tt{u) = e-~ zu ~ l ^\[u + z) . 

In concrete terms, on the standard orthonormal basis, 



Q 



"0 


1 








..." 




' 


1 











1 





V2 





... 




-1 





V2 











V2 





V3 


... 







-V2 





Vs 











Vs 





41 ... 










-Vs 





41 











v 7 ! 


... 



















These matrices satisfy 



QP - PQ = 2i ■ 1 and \{Q 2 + P 2 ) 



H 



where 



H 
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Proof. It only remains to check the matrices for Q and P. We note that 

e^vr(M) = W H (iy)n(u) = e iyu -^ y \(u + iy) , 
and we find by differentiation 

Q7l(u) = U7l(u) + 7t'(m) . 

Taking the coefficient of u n the matrix of Q is found. The matrix for P is found in 
the same way. The choice of the ground state energy c = \% = 1 in the definition of 
H fixes the relation with Q and P correctly. □ 
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The Gaussian representation. Here is another useful representation of the harmonic 
oscillator algebra on a Hilbert space. 

Let H G := L 2 (M,7), where 7 is the standard Gauss measure on R: 

1 _l 2 

7(cfx) := r y(x)dx := e ? x dx . 

V 27T 

Define for zeC the vector by 

: a; i— > e 



Then with 2 6 C is a total set in 7^. (Actually, 2 G zR is already sufficient by 
the uniqueness of the Fourier transform.) Again we find 

(e(z),e(u)) = e zu , (z,ueC). 

Proposition 2.4. There exists a unitary map Una '■ 7~Lh — > 7~£g su °h that for all 
zeC 

Uhg^(z) = e(z) . 

This map sends the vector e n := (0, ■ ■ • , 0, 1, 0, ■ ■ • ) into the n-th Hermite polynomial, 
where these polynomials are given by the generating function 



n=0 

Consequently, this version of the Hermite polynomials satisfies 

1 „ 



h n (x)h m (x) 7(cfe) 



nl 



Proof. The map Tt(z) ^ e(z) extends to a unitary map since the linear spans of the 
ranges of tt and e are dense and both 7r and e are minimal dilations of the positive 
definite kernel (z, u) 1— > e zu . □ 

Let us carry over the relevant operators with this unitary transformation. We find: 

(e isQc ^) (x) = e isx ^(x), (Q&l))(x) = xi/>(x) ; 

(e- itPa i)) (x) = ip(x - 2t) f ^ x ~^A (PgiI>)(x) = ixif)(x) - 2ip'(x) ; 

V 70) / 

H G ^ = (2N G + 1)^ = + 2x— i> + V • 

ox ox 

The Schrddinger representation. Finally we get to the standard Schrodinger repre- 
sentation (|2.1jl . (|2.Hjl and 1)2.6)1 of the harmonic oscillator by dividing away a factor 
\/7(x). Let Us := L 2 (R) and define 

U GS :H g ^H s : {U gs ^){x) := ^W)^) ■ 
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The problem of damping. A damped harmonic oscillator is an evolution (T t ) t >o 
on the real-linear span of a canonical pair (P, Q) that has the form 

T t (Q) = e~ vt (Q cos ut + P sin out) , 

T t (P) = e~ nt (-Q sin cut + P cos cut) , (77 > 0) . (2.16) 

(We apologise for a clash of notation: T t is not related to translations.) This 
spiralling motion in the plane compresses areas by a factor e~ 2rit , so that for t > 
the operators T t {Q) and T t (P) disobey the canonical commutation relation, and T t 
cannot be extended to an automorphism of B{7i). 

Yet this damped oscillatory behaviour occurs in nature, for instance when an atom 
is loosing its energy to its surroundings by emission of light. So it would be worth 
while to make sense of it. There are two basic questions related to this model. 

Question 1. How should T t be extended to £>(7Y)? 

Question 2. Can (T t )t>o be explained as part of a larger whole that evolves by 
*-automorphisms of the form a t (A) = UtAUf 1 , where Uf satisfies a Schrodinger 
equation? 

Spirals and jumps. In Heisenberg's matrix mechanics atoms were supposed to move 
in a mixture of two ways. Most of the time they were thought to rotate according 
to the evolution C/ t as described above, but occasionally they made random jumps 
down the ladder of eigenvalues of the energy operator H. Each time an atom made 
such a jump, it emitted a quantum of light whose (angular) frequency uj was related 
to the size E of the jump by 

E — hijj . 

The probability per unit of time for the atom to jump was given by Fermi's 'Golden 
Rule', formulated in terms of the coupling between the atom and its surroundings, 
and it is proportional to the damping rate r\. 

In the following sections we shall describe this behaviour as a quantum Markov 
process. Both jumps and spirals will be visible in the extension of our T t to the 
atom's full observable algebra. This will be our answer to Question 1, for which we 
shall need the notion of completely positive operators. 

Our answer to Question 2 will be a reconstruction of the atom's surroundings: a 
dilation. There we shall see how the atom can absorb and emit quanta. 

3. Conditional expectations and operations 

We shall now give a sketch of the operational approach to quantum probability 
which was pioneered by Davies, Lewis and Evans ( |Davj . |EvL| ). 

Conditional expectations in finite dimension. In this section we choose for 
definiteness: A := M n , the algebra of all complex n x n matrices, and 

if : A -»• C : A i-> tr (pA) , 

where p is a symmetric n x n matrix with strictly positive eigenvalues and trace 1, 
so that if is faithful. 
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Let A be a symmetric n x n matrix with the spectral decomposition 

aGsp(A) 

The orthogonal projections E a , a G sp(A), form a partition of unity. Measuring 
the observable A means asking all the compatible questions E a at the same time. 
Precisely one of the answers will be 'yes', as stipulated in the interpretation rules. 
If the answer to E a is 'yes', then A is said to take the value a. This happens with 
probability (p(E a ). 

It is natural to define the expectation of A as 

QSsp(A) yoSsp(A) J 

So the state (p not only plays the role of a probability measure, but naturally extends 
to the associated expectation. 

Now let B = B* G A be a second observable with spectral decomposition 

B= 

/3Gsp(B) 

If we first measure B and then A in each trial, in the limit of increasingly many 
trials we obtain a probability measure P on sp(v4) x sp(-B). By the discussion of 
interpretation of quantum probability in Section the probabilities are given by 

P({(a, 0)}) = viFpEaFp) . 

It is then natural to define the conditional probability PL4 = a\B = f3\ as that 
proportion of the trials that have yielded B = (3 which turn out to give A = a later: 

= a\B = 0\:=- P (ft '#» _<f(FpE a F f3 ) 



E Qesp ( 4 )F({(«^)}) <p{Fp) 
The associated conditional expectation is naturally defined as 

E(A\[B = 0\):= £ *nA = a\B = P] = ^j^, 



oesp(A) 



Note that this is a function, / say, of (3. Seen as a quantum random variable this 
conditional expectation is described by the matrix f{B): 

mm) = j( B ) = y, nm= E ,£ ^r F ^ («> 

Note that 

V (E(A\B))= v{FpAFp). 
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Remark. In general we do not have 

<p(E(A\B)) = ip(A) . (3.2) 

The left hand side is the expectation of A after measuring B. The right hand side 
is the expectation of A without any previous operation. The fact that these two 
expectation values can differ is typical for quantum probability. 

Let us give a simple counterexample to (|3.2|) here: Let A := M 2 , choose A G (0, 1), 
and put 



P 





1 - A 



A 



B 



It is readily checked that 

<p{A) = \, <p(E{A\B)) = \ , 
so that the equality ()3.2|) holds if and only if ip is the trace state. 

The conditional expectation given a discrete random variable. Let B denote the (abelian) 
subalgebra of M n generated by 1 and B, a hermitian matrix. In quantum probability 
theory it is accepted practice ( jl'vl-1 ) to call E(A|5) the conditional expectation of A 
given the algebra B, written Pq{A) only if the equality (|3.2|) does hold. The reason 
is that the definition as it stands does not generalise to observables B with continu- 
ous spectrum, or to noncommutative subalgebras B of A. (The value of <p(M(A\B)) 
changes if the possible values of B are not all distinguished while measuring B: if fix 
and /?2 are not distinguished, their eigenspaces group together into a single subspace, 
and in (j3.1j) the projections Fp 1 and Fp 2 are replaced by the projection Fp 1 + Fp 2 .) 
Note that the projections in B are labeled by subsets of sp(-B): 

= { | V c sp(B) } , 

I3£V 

and that B is the linear span of the projections Fp. 

The following is a finite dimensional version of Takesaki's theorem ( |Takj ) on the 
existence of conditional expectations onto von Neumann subalgebras. 

Theorem 3.1. Let B = B* G M n and let B be the *-algebra generated by 1 and B. 
Let ip : M n — > C : A i— > tv(pA) with p strictly positive and tr (p) = 1. Then the 
following are equivalent. 



(3.3) 



(a) There exists a linear map P : M n —>■ B such that 

VA 6 M n V Fe£(0) : V (FAF) = <p(FP(A)F) . 

(b) There exists a linear map P from M n onto B such that 

(i) P maps positive definite matrices to positive definite matrices. 

(ii) P(l) = 1; 

(iii) if o P = ip; 

(iv) P 2 = P. 

(c) Bp = pB. 

If these equivalent conditions hold, then the linear maps P mentioned in (a) and 
(b) are the same. It is called the conditional expectation onto B compatible with (p. 
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Proof, (a) =>- (b): suppose P : M n — * £> is such that ()3.3|) holds. 

Let A > and decompose P(A) as X]/3g sp (B) a /3-^/3 w hh Fp G £ (23). Then apip(Fp) = 
ip(FpP(A)) = ip(F( 3 P(A)F/ 3 ) = ip(FpAFp) > 0. So ap > for all (3 and P(A) > 0. 

Putting A = 1 in (JUBj) we find that for all (3 G sp(P): <p(P/3) = ^(F /3 P(l)F /3 ) = 
V9(F /3 P(1)). Writing P(l) = E/3e S p(i?) e /9-F/J, we see that ep = 1, hence P(l) = 1. 

By putting P = 1 in ()3.3|) . (iii) is obtained. 

Finally, given A, the element P(A) of B is obviously uniquely determined by ()3.3|) . 
But if A G £>, then P(-A) := A clearly satisfies ()3.3|) . It follows that P is an idempotent 
with range B. 

(b) =>- (c): Make a Hilbert space out of A = M n by endowing it with the inner 
product 

(x,y) v :=<p{x*y) . 

We claim that on this Hilbert space P is an orthogonal projection. Since P is idem- 
potent by assumption (b)(iv), it suffices to show that P is a contraction: 

WPiA)^ < \\A\\r . (3.4) 
Given A G M n , define numbers a^eC and bp > by 

P(A)= a ? F ^ P(A*A)= W 

Then from the positivity property (b) (i) it follows that 

V AeC : P((\-l-A)*(\-l-A)) >0. 
This implies that for all j3 G sp(P) and all A G C, 

|A| 2 - {Jap + Xa^) + bp > , 

from which it follows that 

M 2 < bp , that is P{A)*P{A) < P(A*A) . 

Applying <p to the last inequality and using (iii) yields the statement (|3.4|) . So P is 
an orthogonal projection M n — > 23, that is for all A G M n , 

A - P{A) ± v B . 

This means that for all A G M n : 

(p(AB) = ip(P(A)B) and <p(BA) = ip(BP(A)) . 
But then, since B is commutative, 

<p{BA) = V {BP{A)) = <p(P(A)B) = <p(AB) . 

It follows that 

tr (pPA) = tr (pAB) = tr (BpA) , 

and (c) is proved. 

(c) =>- (a): Suppose that Bp = pB. Then for all P G £(£>) and all A G M n , 
ip(FAF) = tr (pPAP) = tr (FpFA) = tr (pP 2 A) = tr (pFA) = <p(FA) . 
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Therefore, defining P(A) by the r.h.s. of (|3.1|) . and putting F = ^2g &v F@ with 
VCsp(B): 

f(FP(A)F) = ^(f7 ] ^ FF ^ = Y,V( F V AF P) 

= v( F pA) = <f (FA) = ip(FAF) . 

(3GV 

□ 



Operations in finite dimension. Let A and B be finite dimensional von Neumann 
algebras, and let A* and B* denote their duals. A linear map T : A — > B defines by 
duality a linear map T* : B* — > .A*. 

The map T* maps states into states if and only if T is positive, that is maps positive 
elements of A to positive elements of B, and is identity preserving. 

The map T is said to be n-positive if T ® id maps positive elements of A <S> M n to 
positive elements of B ® M n : 



04 



n 



> 



T is called completely positive if it is n-positive for all n G N. In that case T* 
states on 6® M n to states on A ® M n . T is called identity preserving if T( 



®id maps 
1^) = Is- 



Definition. An operation T : A — > B is a completely positive identity preserving 
map. Adjoints of operations will also be called operations. 

The idea is that any physical procedure which takes as an input a state on some 
quantum system described by B, and which turns out a state on a quantum system 
described by A must necessarily be of the above kind. Not all operations in the sense 
of the definition can actually be performed, but certainly nothing else is physically 
possible. Indeed any physical operation on a quantum system A should also define a 
physical operation on A <S> 7t, where 1Z stands for some quantum system not affected 
by the operation. The existence of such an 'innocent bystander' outside our quan- 
tum system A should never lead to the prediction by quantum theory of negative 
probabilities. 

The following example shows that complete positivity is strictly stronger than 
positivity. Let 

I : .U> — : 



a b 




a 


c 


c d 


i — > 


b 


d 



Then T(A*A) = T(A)T(A)* > for all A, but 



I id : M 2 ® M 2 ^ M 2 ® M 2 maps 



"l 








l" 




"l 








cf 














to 








1 




















1 








1 








1 













1 



i.e. it maps a one-dimensional projection to a matrix with eigenvalues 1 and —1. 
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Operations on quantum probability spaces. A quantum probability space (A, (p) 
has a canonical representation on a Hilbert space, called the GNS representation after 
Gel'fand, Naimark and Segal. It is the representation of A on TC V , the Kolmogorov 
dilation of the positive definite kernel 

A x A -> C : (A, B) h-> . 

States which are given by density matrices on this space are called normal states on 
A, and the set of all normal states is denoted by A*. 

When we write T : (^4, <p) — * (£>, ip) , we mean that T is a completely positive 
operator A — > B such that T(l^) = and also ip o T = <p. The latter condition, 
which can equivalently be written as 

r*v = (p , 

ensures that T* maps normal states to normal states. This property is only relevant 
for infinite dimensional von Neumann algebras. When speaking of operations between 
quantum probability spaces we shall always imply that the state is preserved. 

Quantum stochastic processes. Let us now consider the category QP whose ob- 
jects are quantum probability spaces and whose morphisms are operations. 

Lemma 3.2 (Schwartz's inequality for completely positive operators). LetT : (A, <p) - 
(B,ip). Then for all A G A, 

T(A*A) > T(A)*T(A) . 

Proof. Let A be represented on TC. By the positivity of T ® id M , 2 we have for all 
A g A, 

(ip®T(A)ip, (T® id) 

Writing this out we obtain 

(T(A*A) - T(A)*T(A))il>) > . 

□ 

Corollary 3.3. IfT : {A, <p) -> (B, ip) then for all A G A 

<p{T{A)*T(A)) < ip(A*A) . 

This inequality states that T is a contraction between the GNS Hilbert spaces of 
(A,<p) and (B,ip). 

Lemma 3.4. T : (A,<p) — > (B,if>) is an isomorphism in the category QP if and only 
if T : A — > B is a * -isomorphism. 

Proof. (Exercise:) Apply Schwartz's inequality to T and to T _1 . □ 

A random variable (cf. Events and random variables, in Section is an injective 
*-homomorphism 

j : {A, if) -> {A, (p) . 



A -1 




A -1 











^©T(A)V) > 
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A quantum stochastic process ( |AFLj ) is a family (Jt)teT °f random variables indexed 
by time T. Here, T is a linearly ordered set such as Z, R, N or R + . If T = R or R + 
we require that for all A e .4. the curve t t— >■ jt(A) is strongly continuous. 

If T is a group, say Z or R, then the process is called stationary provided that 
jt = T t o jo for some representation t i— > T t of T into the automorphisms of (*4, </?). 

Open system interpretation. We are observing a subsystem with observable algebra 
A of a larger environment with algebra .4 that we cannot see. In the Heisenberg 
picture, the smaller algebra is moving inside the larger one. If to < t\ < ■ ■ ■ < t n is a 
sequence of times, and Ex, E 2 , . . . , E n a sequence of events in A, then 

^{E,)]^) ■ ■ ■j tn _ 1 (E n _ 1 )j tn (E n )j tn _ 1 (E n ^ 1 ) ■ ■ .j t2 {E 2 )j tl {E x )) 

is the probability that E\ occurs at time ti, E 2 at time t 2 , ■ ■ ., and E n at time t n . 
Note the double role played here by the time ordering: Unless some of the questions 
jt k {Ek) recur, that is they lie in jt{A) for different values of t, they must be asked in 
the order dictated by the times 

Stochastic process interpretation. In a classical stochastic process (X t )tej the ran- 
dom variable X t is a different one for different times t, so the events concerning X t 
change in time accordingly. If the process is stationary, X t and X s differ by an auto- 
morphism of the underlying probability space. These observations generalise to the 
noncommutative situation. 

Conditional expectations and transition operators. If we are to describe an 
open quantum system such as the damped harmonic oscillator by an internal dynam- 
ics, say T t : A — > A, without reference to its surroundings, we need to be able to 
keep track of an observable A which starts in A at time zero, during its motion away 
from the algebra A at positive times. That is, we need its conditional expectation. 

In view of the discussion of conditional expectations in finite dimensions, we give 
the following general definition. 

Definition. Let j : (A, <p) — > (A, (p) be a random variable. The conditional expecta- 
tion (if it exists) is the unique morphism P : (A, tp) — > (A, (p) for which 

p o j = id A . 

Without proof we state some properties. 

Proposition 3.5. If P : (A,(p) — > (A, <p) is the conditional expectation with respect 
to j : (A, (p) — > (A, (p), then 

Vb^b^V^ : B 1 P(A)B 2 = P{j(B 1 )Aj(B 2 )) . 

In particular 

^Fe £ (A)^AeA ■ (p{j{F)Aj{F)) = <p(FP(A)F) . 
The second line indicates the connection with Theorem Kill 
Markov processes. Let us now apply the above notion to an open quantum system. 
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Two-time-probabilities. Suppose that for all s G T there exists a conditional expec- 
tation P s with respect to j s . Then the probability for F to occur at time s and E at 
time t > s can be written as 

<p(j s (F)j t (E)j s (F)) = <p(FP s (j t (E))F) = <p(FT s>t (E)F) , 

where T s t — P s o j t is an operation on (A, (p), the transition operator from time s to 
time t. 

Multi-time-probabilities. This reduction to the subsystem succeeds for more than two 
time points if there also exist conditional expectations P(-oo,t] onto the algebras 

A-oo,t] : = vN { 3s{A) \ s<t) . 
and moreover the Markov property holds: 

t<s P { ^ t] (js(A)) C j t (A) . (3.5) 

Proposition 3.6. Let (j t : (A,(f) — ► (A, <fi)) t T a Markov process with conditional 
expectations Pt- Then the transition operators form a monoid: 

< s <t < u =>• r S)t r t)U = T SjU . 

In particular, if the process is stationary, then T t := T Q t = T StS+t satisfies 

T s T t = T s+t (s,t>0). 

In the latter case, (T t ) t > is known as the dynamical semigroup induced by the 
stationary Markov process. Conversely, the process (jt)ter is called a Markov dilation 
( |Kiimj ) of the dynamical semigroup (T t ) t6 T- 

The situation is symbolised by the commutative diagram 

{A, (f) {A, cp) (3.6) 

3 P 

^ rp ^ 

{A, tp) {A, <p) 

Our goal is to describe a Markov dilation of the damped harmonic oscillator. 

4. Second quantisation 

A quantum model of n harmonic oscillators is obtained by taking the n-fold tensor 
product of the representation z h- > W(z) := exp(i(Imz)Q— z(Rez)P) of the canonical 
commutation relation (CCR) over C. This turns out to be equivalent to a single 
representation of the CCR over C n . An infinity of harmonic oscillators is obtained by 
replacing C with an infinite dimensional separable Hilbert space /C. It depends on the 
spectrum of the time evolution on K, (discrete or continuous), whether a countable 
infinity of oscillators is obtained or a continuum, that is a quantum field. In our 
dilation of the damped harmonic oscillator we shall need a quantum field. 

As in the case of a single oscillator we have the choice between different concrete 
representations: we may emphasise the field aspect of the construction, like in the 
Gaussian representation of the harmonic oscillator, or the particle aspect of it, like 
in its matrix representation. (The Schrodinger representation on L 2 (M) as in (|2.2|) . 
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()2.6|) has no analogue in infinite dimension, since there exists no Lebesgue measure 
on R°°.) 

The following definition generalises the Weyl relation ()2.14|) over C to that over a 
general complex Hilbert space K. If fC is the L 2 -space of some measure space (X, /i), 
then K, may be considered as the 'quantisation' of X, and the construction below as 
its 'second quantisation'. 

We refer to Mark Fannes' lectures in these volumes. 

The functor T. 

Definition. Let /C be a complex Hilbert space. A representation of the Canonical 
Commutation Relations (CCR) over K, is a map W from K to the unitary operators 
on some Hilbert space Ti such that for all f,ge)C: 

W(f)W(g) = e- am ^W(f + g) , (4.1) 

and t I— > W(tf)ip is continuous for all / E K,, if) £ Ti. The map is called a vacuum 
representation if there is a unit vector Q 6 Ti such that 

(Q,W/(/)Q) = e-^ll 2 . 

A vacuum representation is called cyclic if the linear span of the vectors W(f)Q is 
dense in Ti. 

A cyclic vacuum representation of the CCR over /C can be constructed by a gener- 
alisation of the method used several times in the harmonic oscillator in Sectional Let 
7r be a minimal Kolmogorov decomposition of the positive definite kernel /C x /C — > C 
mapping (/, g) to e^' 9 ', and on the total set of 'coherent vectors' ir(g), define 

W(f)7r(g) := c -</*>-ill/ll a 7r (/ + « / ) , (f,ge)C). 

Then put Q := 7r(0), and all the requirements in the above definition are met. On the 
other hand, given any cyclic vacuum representation of the CCR over JC with vacuum 
vector Q, a Kolmogorov decomposition n' of the above mentioned kernel is obtained 
by putting 

Indeed, 

(ir'(f), n'(g)) = e *<« / " 2+ H ^ • e jIm ^ (O, + y)n> 

= gidl/f + IMI 2 ) . e »Im(/,9> e -|||/-9|| 2 

Thus all cyclic vacuum representations of the CCR over a Hilbert space /C are unitarily 
equivalent. However, this can not be concluded from von Neumann's uniqueness 
theorem, since the latter breaks down for infinite dimesional /C. In this case there are 
indeed many inequivalent (non- vacuum) representations, for instance those associated 
to positive temperatures ( BrRj). 

Since e^ 1 ®^ 2 ' 91 ® 92 ^ = e^ 1 ' 91 ' ■ e^ 2 ' 92 ', a representation of the CCR over a direct sum 
JCi © K.2 of Hilbert spaces is isomorphic to the tensor product of the representations 
of the CCR over /Q and /C 2 . 
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Definition. Let r (/C) denote the linear span of the operators W(f), f G K. in some 
representation of the CCR over /C. Let T(/C) be its strong closure. On the von 
Neumann algebra T(/C) we assume by default the vacuum state 



n2 



Thus r(/C) = (r(/C),<y9^) is a quantum probability space. 

If C is a contraction /Ci — > /C 2 , let r (C) : r (/Ci) — * r (/C 2 ) be given by 

r (C) := e*<« C/ 1|2 -" / » V (C/) • (4.2) 

Proposition 4.1. T/ie operator T (C) has a unique strongly continuous extension to 
an operation T(C) ofT{K). 

Proof. Cf. for instance |Teej . □ 

Remark. Second quantisation is a functor V from the category of Hilbert spaces 
with contractions to the category of quantum probability spaces with operations. 

Fields. From the Weyl relation (|4.1|) it follows that A i— > W(Xf) is a strongly contin- 
uous unitary representation of R. By the spectral theorem there exists a self-adjoint 
operator $(/) on T(/C) such that 

W(Xf) = e lA * (/) . 

The Weyl relations then imply that 

mf),$(g)]=2tim(f,g)-l, 

and in the vacuum state <£>ac the random variable $(/) has normal distribution with 
mean and variance || / 1| 2 . The random variables $(/) and <&(g) are compatible if 
the inner product (/, g) is real, and independent if it is zero. 

In particular, if /C = L 2 (R, 2r]dx) (as we shall need in Sections El and EJ), then by 
putting 

r*(i ;o ,) if(a0 - («) 

\-*(l |t ,o]) if t < 0, 

a stochastic process (B t ) te ^ is defined with compatible normally distributed indepen- 
dent increments having variance 

<pjc((B t -B s ) 2 )= \t-s\. 

Thus B t is a classical Brownian motion. 

Particles. A natural choice for the representation space of the CCR over /C is the 
Fock space J-(K) from Example (d) at the end of Section ^ 



oo 



n=0 
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with 7r(/) given by the exponential vectors. Given a contraction C : /C — > JC let J~{C) 
be the contraction .F(/C) — >■ J-(JC) mapping 7r(/) to n(Cf) for every / G /C. This 
map can be written 

oo 
n=0 

Given an orthogonal projection P on K, let a self-adjoint operator dJ z {P) on !F(K.) 
be defined by 

This operator dJ-'(P) is interpreted as the random variable that counts for how many 
particles the 'question' P is anwered 'yes'. In particular the total number of particles 
N equals 

If /C = L 2 (R, 2r)dx), the Fock space J-{]C) can be written as L 2 (A(R), where 
A(R) is Guichardet's space over R ( |Guij . |Maaj - - see also the lectures by Martin 
Lindsay in these volumes). 

A(R) := { a C R | #(<r) < oo }; 

and n n is the measure on A given by 

^({0}) = 1 , 

H v {da) = (27/)" dsids 2 ■ ■ ■ ds n if a — {si, s 2 , • • • , s n }- 
The coherent vectors are represented as the functions 7r(/) : A(R) — > C given by 

.SG(T 

Indeed, 



W)M9))= / n(f)(a)n(g)(a)^(da) 
Ja(r) 

oo „ 

= ^T(2r/)" / (fg)(s 1 )---(fg)(s n )ds 1 ---ds n 
= E^T / Qg){s 1 ).-.Qg){s n )ds l .--ds n 

°° i / />oo \ n 

= E- 2 W = e(/,s> - 

n=0 n ' V. / 

In this concrete representation the number operator dJ-(P), where P is the mul- 
tiplication in L 2 (R) by Is, counting the number of particles in the region S C R, is 
itself a multiplication operator, multiplying by the number #(cr n S) of quanta in S. 
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This is seen by the following calculation: for all A G R, / G K, and a G A( 
(e^V(/)) (a) = vr ( e ^/) (a) 
= JJ e ^W/( t ) 



tea- 



exp h'A|]ls(*)j 



7r(/)(a) = e lA# ( CTnS ) • 7r(/)(a) 



5. Unitary dilations of spiraling motion 

In preparation for the solution of the physical problem of damping posed at the 
end of Section we now consider embeddings of the spiraling evolution (|2.1fij) into 
a unitary one. Let us describe the spiral by 

C t : C C : z e^+^z , (t > 0). (5.1) 

Theorem 5.1 ((Sz. Nagy, Foias 1953; special case.). Up to unitary equivalence there 
exists a unique Hilbert space K with a unit vector v and a one-parameter group of 
unitary transformations U t on K such that the span of the vectors U t v, t eR is dense 
in /C and 

(v, U t v) = e^+^t , (t>0). 
Proof. Existence: Take JC := L 2 (M, 2r)dx), let U t be the shift to the right, and 



v(x) :- 



if x > 0, 

e fo-tw)t if x < o. 



Then we arrive at Example 11.61 (e) at the end of Section ^ Uniqueness follows from 
Theorem 11.51 □ 

The structure (K.,J,Ut), illustrated in the diagram below, is called a minimal 
unitary dilation of (Ct)t>o- 

e (-»j+iw)t 

C (5.2) 



J* :kt— >(v,k) 



U t 

K -+K 

In practice several — unitarily equivalent — minimal unitary dilations of {C t )t>o can 
be useful. If K, = L 2 (R) and U t is the shift, then we speak of translation dilations of 
(Ct)t>o- They differ only in the shape of v G £ 2 (JR), which must satisfy 

l^)| 2 = 7T ( AeR )- 

(A — UJ) Z + Tf 

Particular solutions are v±(\) :— 1/ (A — u±ir}). Here v + , which occurred in the proof 
of Theorem I5.1[ leads to the incoming translation dilation and t>_ to the outgoing 
translation dilation: 

-(r,+iw)x (a; > 0); 

(x< 0). 



V-(x) 
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The former is more useful for the study of incoming fields and particles, the latter 
for outgoing ones. We shall have occasion to employ both below. The unitary 
equivalence of these two unitary dilations, asserted by Theorem 15.11 is implemented 
by the scattering operator S which in terms of the Fourier transform F can be written 

as 

S:=FM,F-\ 3 (A):=^±^. 

A — uo — ii] 

Apart from these two translation dilations, the interaction dilation (/C, J, Ut), where 
/C = L 2 (— oo, 0] ©C©L 2 [0, oo), J : z i— > 0©,2©0, and U t describes a more complicated 
coupling of C to an incoming and an outgoing channel, is physically more enlightening, 
but too cumbersome to treat here. We refer to jKiiSj for a thorough treatment. 

6. The damped harmonic oscillator 

Equipped with the notions introduced in Sections [TJ El IU and El we are now in a 
position to answer the questions posed at the end of Section El 

We act with the second quantisation functor V of Section 0] on all four corners 
and all four arrows of the dilation diagram (|5.2|) of Section El The corners become 
quantum probability spaces (Section HJ), and the arrows become operations (Section 

r(C) — r(c) (6.i) 



r(j) 



r(j*) 



r(/c) — r(/c) 

The answers to Questions 1 and 2 can now be read off. 

(1) r(C) = BiU), where H = H H = / 2 (N, ±) or equivalently H = H G = L 2 (R, 7) 
is the Hilbert space of the harmonic oscillator in the Heisenberg or in the 
Gaussian representation. The damped time evolution T t : B(TC) — > B{7i) is 
now given by 

T t := r(c t ) = T (e ( -" +iLj)t ) . 

Then 

T t (W(z)) = e^'^-^Wie^+^z) . (6.2) 

By substituting W(t + is) = e - ltp + ts Q an d differentiating with respect to t 
and s respectively, we indeed obtain the equations ()2.16j) . 

(2) The diagram shows how T t is embedded into a larger whole, where the time 
evolution is a one-parameter group of *-automorphisms, that is reversible. 
Here j := T(J) is an injective *-homomorphism, P = T(J*) is a conditional 
expectation. By Theorem 15. II this is the only quasifree dilation, that is in the 
range of the functor T. It is automatically Markov. 

In this Section we shall discuss four aspects of the construction: the stochastic 
behaviour of the oscillator (spirals), its driving field (a quantum Brownian motion), 
the jumps between the levels of the oscillator (a death process), and the outgoing 
quanta (a point process). A complete picture would include the outgoing field and 
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the scattering of incoming particles as well. This can easily be achieved using the 
tools developed here. 

Stochastic behaviour of the oscillator. By the functorial character of T we can 
split T t as 

T t = r(e ( " , ' + ^ )t ) = r(e^*)r(e- r '*) . 

The operator T(e~ tujt ) is the automorphism at studied in Section 01 So let us now 
look at the 'dissipative' part T(e~ vt ). 

Proposition 6.1. For < c < 1 the operator T(c) leaves invariant the abelian 
subalgebras generated by 1 and any of the operators xP — yQ with (x,y) G M. 2 . In 
particular its action on the algebra 

Q ■= { f(Q) | / e L°°(R) } 

is given by 



r(c)(/(Q)) = / e »ci=3) f( c Q + x)dx. (6.3) 

V27r(l - c 2 ) 



Proof. Obviously, T(c) leaves the linear span of j W(Az) | A G K } invariant, and thus 
also its strong closure by Proposition 14.11 Putting f(Q) = e tyQ the r.h.s. of ()6.3|) 
equals 

^ " " " x 2 \ 

e~W^?)e ixy dx ■ e lcyQ = e'^-^W '(icy) , 



y27T(l-C 2 ) J 

which is equal to the l.h.s. by the definition (|4.2j) of T. The theorem follows from 
the strong continuity of T(c). □ 

We recognise the semigroup T t of transition operators restricted to Q as the tran- 
sition operators of a diffusion on K with a drift towards the origin proportional to 
the distance to the origin. 

The driving field. Let us consider the dilation of the semigroup T t for u = 0. We 
take the second quantised incoming translation dilation of Section El and substitute 
it into the diagram (jfi.lj) . Let B t be the Brownian motion given by f)4.3jl . and let Q t 
denote the embedded oscillator $(fj). 

Proposition 6.2. TTie embedded oscillator Q t satisfied the integral equation 

Qt~Qs = -V j Qudu + B t -B s . (6.4) 

J s 

This is the integral version of the stochastic differential equation 

dQt = -vQtdt + dB t . 
So we find an embedded Ornstein-Uhlenbeck process in our Markov dilation. 
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Proof. ( |LeTj ) The following equality between functions in K, holds: 

U t V- - U s v^ = -77 J v u du + l [S)t] , (s < t) . 
Acting with $ on both sides of the equation yields (|6.4j) . 



□ 



Quanta. We now concentrate on another abelian subalgebra of T(C), namely the 
algebra of all diagonal matrices in Heisenberg's matrix mechanics. In terms of the 
operator N denoting the number of excitations of the oscillator, this algebra can be 
written as 

J\f := { f{N) | / e Z°°(N) } ~ 1°°{N) . 
This time we need not put uj = 0. Let (p n denote the state — ► C : / i— > f(n). 

Proposition 6.3. The diagonal algebra M is invariant for T t and 

Vn (T t (s N )) = (1 - e- 2rrt (l -s)) n . 



This is the probability generating function of a pure death process [GrS j with the 
generator 



2n 



Proof. (Sketch. Cf. |Teej for the detailed proof.) We can write s N as a weak integral 
over the operators W(z): 



"0 











1 


-1 











2 


-2 











3 


-3 



ir(l-s) 



e -^M 2 w(z)X(dz) 



(6.5) 



where A denotes the two-dimensional Lebesgue measure on C. This relation can be 
checked by taking matrix elements with respect to coherent vectors. Application of 
T t to both sides of ()6.5j) yields for all u, v 6 C, 



(Tc(u),T t (s N )rc(v)) 



■jfl-e-^O—s)) 



Since this expression is not sensitive to the relative phase of u and v, the operator 
T t (s N ) lies in Af. The statement is proved by comparing the coefficients of (uv ) n on 
both sides. □ 



Emitted quanta. Finally, let us see what happens outside the oscillator while it is 
cascading down its energy spectrum. Since we are interested in outgoing quanta at 
positive times, let us now consider the outgoing translation dilation of (T t )t>o and 
represent T(/C) on the Fock space L 2 (A(M)). We denote the number operator dJ r (P Vt ) 
counting the excitations of the oscillator at time t by N t . 
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From Proposition 16.31 it follows that the diagram (|6.1|) can be restricted to the 
subalgebra ~ l°°(N) 



Af ■ 



■Af 



r(£)^£r(£) 



Here 

j : = r(J) : / » f(N ) ; 

P : = Y{J*) : X ^ ((^Xi>® n >)~ G Z°°(N). 

However, since for different times t and s the functions Vt and t> s are neither parallel 
nor orthogonal, the one-dimensional projections P Vt and P Vs do not commute. And 
since for A, fi G R 

the number operators N t and A^ s do not commute either. So the embedded algebras 
jt(Af) with t e R do not generate an abelian subalgebra of T(/C), as was the case for 
the algebras jt(Q) above. 

For every t G R let us consider the following three number operators. 

N t := dJ r (P Vt ), the number of quanta in the oscillator, 

M t := dT ^M 1[too) j the number of quanta that have not yet left the oscillator, 
K t := dT ^M 1( oo t] j the number of outgoing quanta that have left the oscillator. 

Note that the number M t — N t of incoming quanta is not given by a multiplication 
operator, but the number K t of outgoing quanta is. This is due to the fact that we 
are considering the outgoing translation dilation of T t . Note furthermore that the 
operators M t and K s (s, t G R) all commute. 

For positive times the number M t — N t of incoming quanta has expectation in 
the states of the form $ o P g / 1 (N)) which we consider. So we may expect that 
replacing Nt by M t would lead to an embedded classical Markov chain. 

Proposition 6.4. For t > we have the following commuting diagram involving 
abelian von Neumann algebras. 



Af- 



T t 



■Af 



L°° 



p=r(j*) 



a,^)Sl°°(a,^) 
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Proof. For all z, u G C, t > and s G [0, 1] we have 
(tt(u), P o T(^) o j(s> (z)> = (*(«), P tt(«)> 

oo 

= E 



n=0 

oo 



•/A„(R) 

n=o n - \i=i / 




.7 

oo 

(>xi>{^////: / e _2r,r s 1 i t '°°^ r W 

exp | uz | — e ~ 2r?r 

= exp - e^l - s))) . 
By Proposition 16.21 the latter expression is equal to (7r(u),Tt(s N )ir(z)). □ 

Finally, since K t + M t is equal to the total number of quanta, which is the same 
constant for all times, we conclude that a quantum is emitted at precisely the moment 
that the oscillator makes a downward jump. Moreover, these jumps are made at 
independent exponentially distributed random times. 

These phenomena turn out to be natural consequences of damped harmonic motion 
in a noncommutative description. 
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